Ground state bands of the E(5) and X(5) critical symmetries 
obtained from Davidson potentials through a variational 

procedure 

Dennis Bonatsos* 1 , D. Lenis* 2 , N. Minkov^ 3 , D. Petrellis*, P. P. Raychev^ 4 , 

P. A. Terziev t5 

* Institute of Nuclear Physics, N.C.S.R. "Demokritos" 
GR-15310 Aghia Paraskevi, Attiki, Greece 
t Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences 
72 Tzarigrad Road, BG-1784 Sofia, Bulgaria 

Abstract 

Davidson potentials of the form (3 2 + (3q / ' f3 2 , when used in the original Bohr Hamiltonian 
for 7-independent potentials bridge the U(5) and 0(6) symmetries. Using a variational 
procedure, we determine for each value of angular momentum L the value of (3q at which the 
derivative of the energy ratio Rl = E(L)/E(2) with respect to (3q has a sharp maximum, 
the collection of Rl values at these points forming a band which practically coincides 
with the ground state band of the E(5) model, corresponding to the critical point in the 
shape phase transition from U(5) to 0(6). The same potentials, when used in the Bohr 
Hamiltonian after separating variables as in the X(5) model, bridge the U(5) and SU(3) 
symmetries, the same variational procedure leading to a band which practically coincides 
with the ground state band of the X(5) model, corresponding to the critical point of the 
U(5) to SU(3) shape phase transition. A new derivation of the Holmberg-Lipas formula for 
nuclear energy spectra is obtained by-product. 
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1. Introduction 

The recently introduced E(5) [1] and X(5) [2] models are supposed to describe shape 
phase transitions in atomic nuclei, the former being related to the transition from U(5) 
(vibrational) to 0(6) (7-unstable) nuclei, and the latter corresponding to the transition from 
U(5) to SU(3) (rotational) nuclei. In both cases the original Bohr collective Hamiltonian [3] 
is used, with an infinite well potential in the collective /5-variable. Separation of variables 
is achieved in the E(5) case by assuming that the potential is independent of the collective 
7- variable, while in the X(5) case the potential is assumed to be of the form u((3) +14(7). 
We are going to refer to these two cases as "the E(5) framework" and "the X(5) framework" 
respectively. The selection of an infinite well potential in the /3-variable in both cases is 
justified by the fact that the potential is expected to be flat around the point at which a 
shape phase transition occurs. Experimental evidence for the occurence of the E(5) and 
X(5) symmetries in some appropriate nuclei is growing ([4, 5] and [6, 7] respectively). 

In the present work we examine if the choice of the infinite well potential is the optimum 
one for the description of shape phase transitions. For this purpose, we need one-parameter 
potentials which can span the U(5)-0(6) region in the E(5) framework, as well as the U(5)- 
SU(3) region in the X(5) framework. It turns out that the exactly soluble [8, 9] Davidson 
potentials [10] 

= /3 2 + ||, (1) 



where /3q is the position of the minimum of the potential, do possess this property. Taking 
into account the fact that various physical quantities should change most rapidly at the 
point of the shape phase transition [11], we locate for each value of the angular momentum 
L the value of (3q for which the rate of change of the ratio Rl = E(L)/E(2), a widely used 
measure of nuclear collectivity, is maximized. It turns out that the collection of Rl ratios 
formed in this way in the case of a potential independent of the 7-variable correspond to 
the E(5) model, while in the case of the u(/3) +u(j) potential lead to the X(5) model, thus 
proving that the choice of the infinite well potential made in Refs. [1, 2] is the optimum 
one. The variational procedure used here is analogous to the one used in the framework 
of the Variable Moment of Inertia (VMI) model [12], where the energy is minimized with 
respect to the (angular momentum dependent) moment of inertia for each value of the 
angular momentum L separately. 

In Section 2 the E(5) case is considered, while the X(5) case is examined in Section 3, 
in which a new derivation of the Holmberg-Lipas formula [13] for nuclear energy spectra 
is obtained as a by-product. Finally, Section 4 contains a discussion of the present results 
and plans for further work. 

2. Davidson potentials in the E(5) framework 

The original Bohr Hamiltonian [3] is 
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where (3 and 7 are the usual collective coordinates describing the shape of the nuclear 
surface, Qk (k — 1, 2, 3) are the components of angular momentum, and B is the mass 
parameter. 

Assuming that the potential depends only on the variable (3, i.e. V(/3, 7) = U((3), one 
can proceed to separation of variables in the standard way [3, 14], using the wavefunc- 
tion \I/(/3,7,#i) = /( / 9)$(7,^), where 6i (i = 1,2,3) are the Euler angles describing the 
orientation of the deformed nucleus in space. 

In the equation involving the angles, the eigenvalues of the second order Casimir operator 
of SO(5) occur, having the form A = r(r + 3), where r = 0, 1, 2, ... is the quantum number 
characterizing the irreducible representations (irreps) of SO (5), called the "seniority" [15]. 
This equation has been solved by Bes [16]. 

The "radial" equation can be simplified by introducing [1] reduced energies e = ^-E 
and reduced potentials u = jpU, leading to 



1 d nA d t(t + 3) . 



f{(3)=ef{(5). 



(3) 



When plugging the Davidson potentials of Eq. (1) in the above equation, the Pq/P 2 
term is combined with the r(r + 3)//3 2 term appearing there and the equation is solved 
exactly [8, 9], the eigenfunctions being Laguerre polynomials of the form 
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where T(n) stands for the T- function, while p is determined by [8] 

P (p + 3) = t(t + 3) + $, 

leading to 



(4) 



P=- 2 + 
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The energy eigenvalues are then [8, 9] (in Tiu = 1 units) 



E n>T = 2n + p+- = 2n + l + 



1/2 



(7) 



For (3q = the original solution of Bohr [3], which corresponds to a 5-dimensional (5-D) 
harmonic oscillator characterized by the symmetry U(5) D SO(5) D SO(3) D SO(2) [17], 
is obtained. The values of angular momentum L contained in each irrep of SO (5) (i.e. for 
each value of r) are given by the algorithm [18] r = 3z/ A + A, where v A — 0, 1, . . . is the 
missing quantum number in the reduction SO(5) D SO(3), and L — A, A + 1, . . . , 2A — 2, 2A 
(with 2 A — 1 missing). 
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The levels of the ground state band are characterized by L = 2r and n — 0. Then the 
energy levels of the ground state band are given by 



E 0tL = l + - [(L + 3) 2 + 4/3 4 



1/2 



(8) 



while the excitation energies of the levels of the ground state band relative to the ground 
state are 
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For u(f3) being a 5-D infinite well 

«(/?) = 
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one obtains the E(5) model of Iachello [1] in which the eigenf unctions are Bessel functions 
Jt+3/2(z) (with z = j3k, k = y/e), while the spectrum is determined by the zeros of the 
Bessel functions 
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where x^ tT is the £-th zero of the Bessel function J T+3 / 2 (z). The spectra of the E(5) and 
Davidson cases become directly comparable by establishing the formal correspondence n = 

It is instructive to consider the ratios 



Rl = 
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where the notation E n ^ is used. 

For (3 — it is clear that the original vibrational model of Bohr [3] (with R 4 = 2) is 
obtained, while for large (3q the 0(6) limit of the Interacting Boson Model (IBM) [18] (with 
Ra = 2.5) is approached [8]. The latter fact can be seen in Table 1, where the Rl ratios 
for two different values of the parameter (3 are shown, together with the 0(6) predictions 
(which correspond to E(L) = AL(L + 6), with A constant [19]). It is clear that the 0(6) 
limit is approached as /3q is increased, the agreement being already quite good at f3o — 5. 

It is useful to consider the ratios Rl, defined above, as a function of (3 . As seen in Fig. 
1, where the ratios R4, R12 and R20 are shown, these ratios increase with /3q, the increase 
becoming very steep at some value (3o :ma x of (3q, where the first derivative reaches a 

maximum value, while the second derivative ^3r- vanishes. Numerical results for (3o max 
are shown in Table 2, together with the values of Rl occuring at these points, which are 
compared to the Rl ratios occuring in the ground state band of the E(5) model [1]. Very 
close agreement of the values determined by the procedure described above with the E(5) 
values is observed in Table 2, as well as in Fig. 2, where these ratios are also shown, together 
with the corresponding ratios of the U(5) and 0(6) limits. 
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The work performed in this section is reminiscent of a variational procedure. Wishing 
to determine the critical point in the shape phase transition from U(5) to 0(6), one chooses 
a potential (the Davidson potential) with a free parameter (/3 ), which helps in covering 
the whole range of interest. Indeed, for @ = the U(5) picture is obtained, while large 
values of flo lead to the 0(6) limit. One then needs a physical quantity which can serve as 
a "measure" of collectivity. For this purpose one considers the ratios Rl, encouraged by 
the fact that these ratios are well-known indicators of collectivity in nuclear structure [20]. 
Since at the critical point (if any) one expects the collectivity to change very rapidly, one 
looks, for each Rl ratio separately, for the value of the parameter at which the change of 
Rl is maximum. Indeed, the first derivative of the ratio Rl with respect to the parameter 
(5q exhibits a sharp maximum, which is then a good candidate for being the critical point 
for this particular value of the angular momentum L. The Rl values at the critical points 
corresponding to each value of L form a collection, which should correspond to the behaviour 
of the ground state band of a nucleus at the critical point. The infinite well potential used 
in E(5) succeeds in reproducing all these "critical" Rl ratios in the ground state band for all 
values of the angular momentum L, without using any free parameter. It is therefore proved 
that the infinite well potential is indeed the optimum choice for describing the ground state 
bands of nuclei at the critical point of the U(5) to 0(6) shape phase transition. 

In other words, starting from the Davidson potentials and using a variational procedure, 
according to which the rate of change of the Rl ratios as a function of the parameter f3 is 
maximized for each value of the angular momentum L separately, one forms the collection 
of critical values of R L which corresponds to the ground state band of the E(5) model, 
which is supposed to describe nuclei at the critical point. 

Variational procedures in which each value of the angular momentum L is treated sepa- 
rately are not unheard of in nuclear physics. An example is given by the Variable Moment 
of Inertia (VMI) model [12], in which the energy of the nucleus is minimized with respect 
to the (angular momentum dependent) moment of inertia for each value of the angular 
momentum separately. From the cubic equation obtained from this condition, the moment 
of inertia is uniquely determined (as a function of angular momentum) in each case. The 
collection of energy levels occuring by using in the energy formula the appropriate value of 
the moment of inertia for each value of the angular momentum L forms the ground state 
band of the nucleus. 

Some comparison of the variational procedure used here with the standard Ritz varia- 
tional method used in quantum mechanics ([21], for example) is in place. In the (simplest 
version of the) Ritz variational method a trial wave function containing a parameter is 
chosen and subsequently the energy is minimized with respect to this parameter, thus de- 
termining the parameter value and, after the relevant substitution, the energy value. In 
the present case a trial potential containing a parameter is chosen and subsequently the 
rate of change of the physical quantity (here the rate of change of the energy ratios) is 
maximized with respect to this parameter, thus determining the parameter value and, after 
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the relevant calculation, the value of the physical quantity (here the energy ratios). The 
main similarity between the two methods is the use of a parameter-dependent trial wave 
function/trial potential respectively. The main difference between the two methods is that 
in the former the relevant physical quantity (the energy) is minimized with respect to the 
parameter, while in the latter the rate of change of the physical quantity (the energy ratios) 
is maximized with respect to the parameter. 

3. Davidson potentials in the X(5) framework 

Starting again from the original Bohr Hamiltonian of Eq. (2), one seeks solutions of 
the relevant Schrodinger equation having the form 7, 0j) = 0^(/3, K (6i), where 
9i (i = 1, 2, 3) are the Euler angles, £>(#«) denote Wigner functions of them, L are the 
eigenvalues of angular momentum, while M and K are the eigenvalues of the projections of 
angular momentum on the laboratory-fixed z-axis and the body-fixed z'-axis respectively. 

As pointed out in Ref. [2], in the case in which the potential has a minimum around 
7 = one can write the last term of Eq. (2) in the form 



Q 
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(13) 



Using this result in the Schrodinger equation corresponding to the Hamiltonian of Eq. 
(2), introducing reduced energies e = 2BE/H 2 and reduced potentials u = 2BV/H 2 , and 
assuming that the reduced potential can be separated into two terms, one depending on j3 
and the other depending on 7, i.e. u(/3, 7) = u{f3) + ^(7), the Schrodinger equation can be 
separated into two equations [2], the "radial" one being 
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When plugging the Davidson potentials of Eq. (1) in this equation, the (3q/(3 2 term of 
the potential is combined with the L(L + l)/3/5 2 term appearing there and the equation is 
solved exactly, the eigenfunctions being Laguerre polynomials of the form 
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where a is given by 
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The energy eigenvalues are then (in Tiuj = 1 units) 
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The levels of the ground state band are characterized by n = 0. Then the excitation 
energies relative to the ground state are given by 



F - \ l - 

- C/ 0,Z/,exc g 

which can easily be put into the form 
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which is the same as the Holmberg-Lipas formula [13] 



1/2 
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(19) 
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It is clear that the Holmberg-Lipas formula gives rotational spectra for small values of 
bn, at which one can keep only the first L-dependent term in the Taylor expansion of the 
square root appearing in Eq. (20), leading to energies proportional to L(L + 1). From 
Eq. (21) it is then clear that rotational spectra are expected for large values of f3 Q . This 
can be seen in Table 3, where the Rl ratios occuring for two different values of @ are 
shown, together with the predictions of the SU(3) limit of IBM, which correspond to the 
pure rotator with E(L) = AL(L + 1), where A constant [18]. The agreement to the SU(3) 
results is quite good already at (5q — 5. On the other hand, the case f3 = corresponds to 
an exactly soluble model with R4 = 2.646, which has been called the X(5)-/3 2 model [22]. 

It is worth remarking at this point that the Holmberg-Lipas formula can be derived 
[19] by assuming that the moment of inertia I in the energy expression of the rigid rotator 
(E(L) = L(L + 1)/2J) is a function of the excitation energy, i.e. I — a + f3E(L), where a 
and (3 are constants, the latter being proportional to b H and acquiring positive values. It is 
therefore clear that the Holmberg-Lipas formula, as well as the spectrum of the Davidson 
potentials derived in this section, have built-in the concept of the Variable Moment of Inertia 
(VMI) model [12], according to which the moment of inertia is an increasing function of 
the angular momentum. 

For u(f3) being a 5-D infinite well potential (see Eq. (10) one obtains the X(5) model 
of Iachello [2], in which the eigenfunctions are Bessel functions J v {k s ^f3) with 



' L(L + 1) | 9 s 



1/2 
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while the spectrum is determined by the zeros of the Bessel functions, the relevant eigen- 
values being 

e/3 ;s ,L = (k s ,L) 2 , k s ,L = (23) 

Pw 

where x s ^l is the s-th zero of the Bessel function Ju(k s ^P)- The spectra of the X(5) 
and Davidson cases become directly comparable by establishing the formal correspondence 
n — s — 1. 

It is useful to consider the ratios Rl, defined in the previous section, as a function of 
(5q. As seen in Fig. 3, these ratios again increase with f3 , the increase becoming very steep 
at some value (3^ max of fio, where the first derivative reaches a maximum value, while 
the second derivative ^Mr- vanishes. Numerical results for (3o max are shown in Table 4, 
together with the values of Rl occuring at these points, which are compared to the Rl 
ratios occuring in the ground state band of the X(5) model [2]. Very close agreement of the 
values determined by the procedure described above with the X(5) values is observed. 

The work performed here is reminiscent of a variational procedure, as in the previous 
section. Wishing to determine the critical point in the shape phase transition from U(5) to 
SU(3), one chooses a potential (the Davidson potential) with a free parameter (/3q), which 
serves in spanning the range of interest. For large values of /3q the SU(3) limit is obtained, 
while for (3 = the X(5)-/3 2 picture is obtained [22], which is not the U(5) limit, but it is 
located between U(5) and X(5), on the way from U(5) to SU(3). Thus the region of interest 
around X(5) is covered from X(5)-f3 2 to SU(3). Then the values of (3 at which the first 
derivative dRi/d^o exhibits a sharp maximum are determined for each value of the angular 
momentum L separately, the collection of Rl ratios at these values of (3q forming a band, 
which turns out to be in very good agreement with the ground state band of X(5), the model 
supposed to be appropriate for describing nuclei at the critical point in the transition from 
U(5) to SU(3), thus indicating that the choice of the infinite well potential used in the X(5) 
model is the optimum one. The results are depicted in Fig. 4, where in addition to the 
bands provided by the variational procedure and the X(5) model, the bands corrsponding 
to the U(5), X(5)-/5 2 , and SU(3) cases are shown. 

4. Discussion 

The main results and conclusions obtained in the present work are listed here: 

1) A variational procedure for determining the values of physical quantities at the point 
of shape phase transitions in nuclei has been suggested. Using one-parameter potentials 
spanning the region between the two limiting symmetries of interest, the parameter values 
at which the rate of change of the physical quantity becomes maximum are determined 
for each value of the angular momentum separately and the corresponding values of the 
physical quantity at these parameter values are calculated. The values of the physical 
quantity collected in this way represent its behaviour at the critical point. 

2) The method has been applied in the shape phase transition from U(5) to 0(6), using 
one-parameter Davidson potentials [10] and considering the energy ratios Rl = E(L)/E(2) 
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within the ground state band as the relevant physical quantity, leading to a band which 
practically coincides with the ground state band of the E(5) model [1]. It has also been 
applied in the same way in the shape phase transition from U(5) to SU(3), leading to a 
band which practically coincides with the ground state band of the X(5) model [2]. 

3) It should be emphasized that the application of the method was possible because the 
Davidson potentials correctly reproduce the U(5) and 0(6) symmetries in the former case 
(for small and large parameter values respectively), as well as the relevant X(5)-/3 2 [22] and 
SU(3) symmetries in the latter case (for small and large parameter values respectively). 

4) As a by-product, a derivation of the Holmberg-Lipas formula [13] has been achieved 
using Davidson potentials in the X(5) framework. 

It is clearly of interest to apply the variational procedure introduced here to physical 
quantities other than the energy ratios in the ground state band. Energy ratios involving 
levels of excited bands, ratios of B(E2) transition rates (both intraband and interband), and 
ratios of quadrupole moments are obvious choices. Work in these directions is in progress, 
using the Davidson potentials, since they possess the appropriate limiting behaviour for 
small and large parameter values. However, any other potential/Hamiltonian bridging the 
relevant pairs of symmetries (U(5)-0(6) and U(5)-SU(3)) should be equally appropriate. 
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Figure captions 

Fig. 1 The R L ratios (defined in Eq. (12)) for L — 4, 12, 20 and their derivatives 
dR,L/d[3o vs. the parameter (3 , calculated using Davidson potentials (Eq. (1)) in the E(5) 
framework. The Rl curves also demonstrate the evolution from the U(5) symmetry (on 
the left) to the 0(6) limit (on the right). See section 2 for further details. 

Fig. 2 Values of the ratio Rl (defined in Eq. (12)) obtained through the variational 
procedure (labeled by "var") using Davidson potentials in the E(5) framework, compared 
to the values provided by the U(5), 0(6), and E(5) models. See section 2 for further details. 

Fig. 3 The R L ratios (defined in Eq. (12)) for L — 4, 12, 20 and their derivatives 
dR L /df3 vs. the parameter (3 , calculated using Davidson potentials (Eq. (1)) in the X(5) 
framework. The Rl curves also demonstrate the evolution from the X(5)-/5 2 symmetry (on 
the left) to the SU(3) limit (on the right). See section 3 for further details. 

Fig. 4 Values of the ratio Rl (defined in Eq. (12)) obtained through the variational 
procedure (labeled by "var") using Davidson potentials in the X(5) framework, compared 
to the values provided by the U(5), SU(3), X(5), and X(5)-/3 2 models. See section 3 for 
further details. 
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Table 1: R L ratios (defined in Eq. (12)) for the ground state band of the Davidson potentials 
in the E(5) framework (Eq. (8)) for different values of the parameter (3 , compared to the 
0(6) exact results. 



L 


Rl 


Rl 


Rl 




A) = 5. 


A) = io. 


0(6) 


4 


2.494 


2.500 


2.500 


6 


4.475 


4.498 


4.500 


8 


6.935 


6.996 


7.000 


10 


9.861 


9.991 


10.000 


12 


13.242 


13.483 


13.500 


14 


17.064 


17.471 


17.500 


16 


21.312 


21.954 


22.000 


18 


25.969 


26.930 


27.000 


20 


31.020 


32.398 


32.500 



Table 2: Parameter values Po,max where the first derivative of the energy ratios Rl (de- 
fined in Eq. (12)) in the E(5) framework has a maximum, while the second derivative 
vanishes, together with the R L ratios obtained at these values (labeled by "var") and the 
corresponding ratios of the E(5) model, for several values of the angular momentum L. 



L 


Po,max 


Rl 


Rl 






var 


E(5) 


4 


1.421 


2.185 


2.199 


6 


1.522 


3.549 


3.590 


8 


1.609 


5.086 


5.169 


10 


1.687 


6.793 


6.934 


12 


1.759 


8.667 


8.881 


14 


1.825 


10.705 


11.009 


16 


1.888 


12.906 


13.316 


18 


1.947 


15.269 


15.799 


20 


2.004 


17.793 


18.459 
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Table 3: Rl ratios (defined in Eq. (12)) for the ground state band of the Davidson potentials 
in the X(5) framework (Eq. (17)) for different values of the parameter /3 , compared to the 
SU(3) exact results. 



L 


Rl 


Rl 


Rl 




A) = 5. 


A) = io. 


SU(3) 


4 


3.327 


3.333 


3.333 


6 


6.967 


6.998 


7.000 


8 


11.897 


11.993 


12.000 


10 


18.087 


18.317 


18.333 


12 


25.503 


25.968 


26.000 


14 


34.102 


34.941 


35.000 


16 


43.839 


45.233 


45.333 


18 


54.665 


56.841 


57.000 


20 


66.530 


69.760 


70.000 



Table 4: Parameter values Po,max where the first derivative of the energy ratios Rl (de- 
fined in Eq. (12)) in the X(5) framework has a maximum, while the second derivative 
vanishes, together with the R L ratios obtained at these values (labeled by "var") and the 
corresponding ratios of the X(5) model, for several values of the angular momentum L. 



L 


Po,max 


Rl 


Rl 






var 


X(5) 


4 


1.334 


2.901 


2.904 


6 


1.445 


5.419 


5.430 


8 


1.543 


8.454 


8.483 


10 


1.631 


11.964 


12.027 


12 


1.711 


15.926 


16.041 


14 


1.785 


20.330 


20.514 


16 


1.855 


25.170 


25.437 


18 


1.922 


30.442 


30.804 


20 


1.985 


36.146 


36.611 
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